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We search the COSMOS survey for pairs of galaxies consistent with the gravitational lensing 
signature of a cosmic string. The COSMOS survey imaged 1.64 square degrees using the Advanced 
Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST). Our technique includes 
estimates of the efficiency for finding the lensed galaxy pair. We find no evidence for cosmic strings 
with a mass per unit length of Gji/c 2 < 3.0 x 10 -7 out to redshifts greater than 0.6 and set 95% 



upper limits. This corresponds to a global 95% upper limit of Q 



strings 



< 0.0028. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION: 



Cosmic strings are linear topological defects that arise 
naturally during symmetry-breaking phase transitions in 
the early universe p], Q . They have also been proposed 
in string theory models of inflation, occuring just after 
the GUT scale transition [3]. Detailed simulations of 
the dynamics and interactions of cosmic strings predict 
a modern day stochastic network of strings, observable 
through a variety of astrophysical phenomena [4|, [5[ . The 
dimensionless scale of observational interest is G/i/c 2 < 
10 -6 where \i is the string energy-density [fj. 

Despite considerable interest within the theory com- 
munity and multiple proposed production mechanisms, 
only a few observations bear on the subject, including 
cosmic microwave background (CMB) [7Hlll|, gravita- 
tional waves jl2l - [l4) , and gravitational lensing [15|, [l6| . 
The CMB power spectrum shows that cosmic strings are 
not the dominant factor in large-scale structure forma- 
tion, but that they may contribute up to about 10% of 
the observed structure, enough to possibly detect with 
the Planck data [13, [H| or other experiments with reso- 
lution on small angular scales . 

Gravitational lensing of background galaxies by a cos- 
mic string is expected to produce a pair of images sep- 
arated by an angle, A6 = 5 sm(f3)Di s /D os where S — 
87rGfi/c 2 is the deficit angle, Di s is the distance between 
the lensing string (I) and the background source (s), D os 
is the distance between the observer (o) and the back- 
ground source, and f3 is the tilt of the string toward the 
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observer. 

In our previous search for gravitational lensing by cos- 
mic strings in the GOODS survey we concluded with 
95% confidence that Gfi/c 2 < 3.0 x 10 -7 out to red- 
shifts greater than 0.5 and that £l s trings < 0.02 [15]. Our 
aim in this paper is to use the same technique to analyze 
the wider survey carried out by the COSMOS team with 
the Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS). This survey has a field-of-view that 
is 24 times larger than the GOODS survey and the limits 
we find are nearly 10 times lower than our previous lim- 
its. Because strings are linear objects and backgrounds 
scale with the field-of-view, we have found it necessary 
to perform the search in four angular bins to improve the 
signal-to-noise ratio. 

In Section II we describe our data selection and present 
the correlation analysis used to search for cosmic strings 
in the COSMOS survey. We discuss the simulations 
needed to estimate the signal rates and detection effi- 
ciencies in Section III. These estimates are then used 
in Section IV to determine limits on individual cosmic 
strings as a function of mass and redshift as well as the 
global limit on the density of cosmic strings. Finally, we 
summarize our results in Sec. V. 



II. DATA SAMPLE: 

The COSMOS field is 1.64 degrees 2 centered on 
RA=10:00:28.6 and DEC=+02:12:21.0. Images were 
taken with the ACS aboard HST between July 2003 and 
June 2005 [20|, Hl|. We analyze the publicly available 
COSMOS Version 1.3 data in the F814W (I-band) filter 
which consists of 81 drizzled tiles with a resolution of 
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0.05 "/pixel. We apply a fiducial cut and use the central 
1.57 degrees 2 of the survey. 



A. Source identification 

We use SExtractor version 2.5.0 (Source Extractor) to 
identify sources in the COSMOS survey [22|,l23j following 
the Hot procedure outlined in [24|. For shape-sensitive 
analyses like weak lensing, Leauthaud et al. advocate a 
Hot- Cold method applied to un-rotated, undrizzled im- 
ages. We validate our Hot catalogs against the public 
release of the well-understood Leauthaud et al. cata- 
log containing 1.2 million objects which we refer to as 
the LC2 catalog described in [241 ]. We do not perform 
a catalog- level search using the LC2 catalog directly due 
to the difficulties in estimating the efficiency of finding 
lensed galaxies without access to the original COSMOS 
fits images and processing pipeline. 

Several minor modifications need to be applied to 
Leauthaud's Hot parameters [24] to account for the differ- 
ence in resolution for the publicly available tiles. We set 
the PIXEL_SCALE to 0.5 arcsec, use a gaussian filter width 
of 2.5 pixels, and set the DETECT.MINAREA to 9 contigu- 
ous pixels above threshold. We also scale the apertures 
to 12 pixels, PHOT.APERTURES and PHOT.AUTOAPERS. Fi- 
nally, to reduce the deblending slightly, we set the 
DEBLEND.NTHRESH to 32 and the DEBLEND J4INC0NT to 
0.1. The resulting catalog contains 812,463 objects with 
magnitudes brighter than 26.5. 



B. Resolved galaxy selection 

We select the resolved galaxies from this sample using 
the correlation of the peak surface brightness, MU_MAX, 
with the objects magnitude, MAG_AUT0. Fig.[T]shows three 
regions of interest in the MU_MAX vs MAG_AUT0 plane: re- 
solved galaxies, point sources including stars, and spuri- 
ous detections where the objects are too small to be con- 
sistent with the point spread function. There are 761,370 
resolved galaxies, 36,823 point sources, and 14,270 spu- 
rious detections. The LC2 catalog contains fewer spu- 
rious detections due to a prior cleaning procedure that 
included merging of small objects with nearby objects 
from the Cold catalog and removal of objects in regions 
of elevated noise that occur on the borders of the un- 
rotated tiles as well as near bright stars. To reproduce 
the additional "by hand" cleaning, we correlate our re- 
solved galaxies with the LC2 catalog. 

We find that 96.4% of the resolved galaxies with mag- 
nitudes brighter than 23 rd magnitude are also found in 
the LC2 catalog. Our magnitudes generally agree, but 
have a tail down to smaller values which we attribute to 
over deblending of larger objects. The 3.6% of events 
that don't have a counterpart in the LC2 catalog are all 
found in the diffraction trails of very bright stars or in 
a few cases near the edge of the survey consistent with 
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FIG. 1: Resolved galaxies (black points) are a distinct pop- 
ulation. Also shown are point sources (dark gray) including 
stars, and spurious objects (light gray) that are too small to 
be consistent with the PSF. 



the more aggressive cleaning. The situation is under- 
standably slightly worse for dimmer objects; 87.7% of 
the galaxies at 25 th magnitude have a counterpart in the 
LC2 catalog. The magnitudes tend to agree well. The 
12.3% of galaxies that don't have a counterpart in the 
LC2 catalog are again in fiducial regions that were re- 
moved by hand. For the rest of the galaxies, most of 
them are in regions with elevated noise. Very few appear 
to be legitimate detections. 

To reproduce the cleaning as much as possible, we re- 
move resolved sources in our catalog that are not included 
in the LC2 catalog. This removes about 11.6% of our re- 
solved galaxies. Any inefficiency that comes about from 
this requirement is included in our efficiency estimate. 
More importantly, though, this requirement protects us 
from overestimating the cosmic string lensing rate due to 
spurious detections. 

We post-process our resolved galaxy catalog to identify 
the pixels in the image associated with each galaxy. First 
we define a small but encompassing search region about 
each galaxy centroid. This region is chosen to be three 
times the size of the galaxy reported by the catalog co- 
ordinates (XMIN_ IMAGE, YMIN.IMAGE) and (XMAX_ IMAGE, 
YMAX_ IMAGE). Next we determine the local background 
characteristics by fitting a gaussian to the small ampli- 
tude peak in a histogram of the pixel intensities. Finally, 
we find a bright pixel near the galaxy centroid and it- 
er at ively aggregate neighboring pixels that are la above 
the mean background. This process occasionally merges 
neighboring galaxies. In the event that a cluster of pixels 
reaches the edge of the search region or that two galaxies 
merge, we raise the neighbor threshold to 2a above the 
mean background and repeat the process. We continue 
to raise the threshold until each galaxy is completely con- 
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tained within the search region and does not contain the 
centroid from any other galaxy in the catalog. The aim 
of this procedure is to retain as much unbiased shape in- 
formation as possible. For 0.15% of the dimmest sources 
the threshold is raised so high that there are no pixels 
left in the cluster and we remove these galaxies from the 
sample. After selecting resolved galaxies also identified 
in the LC2 catalog and identifying the galaxy pixels, the 
resulting catalogs contain 662,765 resolved galaxies. 



C. Matched galaxy pair selection 

The morphological similarity between each pair of 
galaxies is characterized by the correlation and cross- 
correlation of the two galaxy images [15j. We first align 
the centroids and then calculate the correlation (CORK) 
and the cross-correlation (XCORR) of the pixel intensi- 
ties. 



CORR = 



XCORR 



(1) 



(2) 



where I(xi,yi) is the intensity of each pixel in a galaxy 
and the subscript 1 or 2 refers to the galaxies being cor- 
related. Galaxies with CORR near a value of zero have 
very similar magnitudes and galaxies with XCORR near 
a value of one are similar in surface brightness and shape. 
We define matched galaxy pairs as those within the el- 
lipse defined by 



V(2 * CORR) 2 + (1 - XCORR) 2 < 0.29 (3) 

This cut was optimized on simulated lensing events and is 
slightly looser than the one used in our previous GOODS 
search. 

In this analysis, we consider pairs of galaxies with 
opening angles, AO < 15 " . There are 96,413 matched 
pairs out of 7,081,011 total pairings with A6 < 15 " in 
the COSMOS survey. The selected matched pairs are 
consistent with the null hypothesis. 



D. Pairs distribution 

The binned distribution of matched galaxy pairs is 
shown in Fig. [2j The distribution is divided into four 
overlapping angular bins: -5°:50°, 40°:95°, 85°:140°, and 
130°: 185°. The angular bins correspond to a range of 
cosmic string angles on the sky. The hypothetical string 
is presumed to pass between the two galaxies in a pair at 
an angle perpendicular to the line connecting the galaxy 
centroids. These bins are needed to reduce the back- 
ground as the survey gets increasingly large. Because 
cosmic strings are linear objects, the signal scales with 
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FIG. 2: Pairs of galaxies in four angular bins (points) are 
compared to background (solid line). The top panel shows 
one example of a simulated string (grey). The upper simu- 
lated string is the total number of pairs expected from the 
simulation with string length of 1.19°, redshift of 0.5, and 
5 sin p of 4 /; . The lower simulated string includes measure- 
ment inefficiencies. 
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the width of the survey whereas the background scales 
with the area of the survey. We do not want to restrict 
our search to perfectly straight strings so we keep the 
number of bins to a minimum. Overlapping bins are 
chosen to assure us that we are searching efficiently. 

The background shape is characterized by the dis- 
tribution of all pairs of galaxies regardless of size and 
shape. Because strings with masses large enough to cre- 
ate opening angles greater than 1" have been ruled out, 
we normalize the background distribution to the number 
of measured matched pairs between 7" and 15 " . This 
gives us a reliable estimate of the background at smaller 
opening angles. From the background, we observe that 
SExtractor merges galaxies with opening angles below 
0.4". 

In our signal region, between 0.4" and 7", there are 
24 data bins and the x 2 /24 — dof of the matched pairs 
to the background is 0.99, 1.1, 1.3, and 1.4 for the four 
angular bins respectively. The p- values range from 15% 
to 48%. Based on the scaled background distribution, we 
report no evidence for an excess of pairs at small open- 
ing angles. Although there is no statistically significant 
excess in our signal region, we note that there is a small 
excess of signal below 1 " and that the first data point 
in the signal region is systematically high. This excess is 
consistent with the fact that the dimmest galaxies tend 
to pass the correlation cuts more easily than larger galax- 
ies. To investigate this region more carefully, we surveyed 
matched pairs with opening angles less than 1 " by hand. 
We find that matched pairs tend to be discovered in re- 
gions of high galaxy density and that these regions tend 
to have other lensing candidates at the fairly high rate of 
0.5/1 ' . A full statistical analysis of these pairs shows no 
evidence of lensing and gives us a calibration curve for 
the background in high density regions. 



III. COSMIC STRING SIMULATIONS: 
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FIG. 3: The simulated redshift distribution for galaxies in a 
typical COSMOS tile. 



of the result from this simulation is shown in the top 
panel of Fig. [2j 

Fig. 3] summarizes the efficiencies as a function of the 
opening angle between the galaxies and the string red- 
shift. These efficiencies were estimated by embedding 
galaxies as realistically as possible into COSMOS tiles 
and then processing the modified tiles back through our 
analysis chain to see how many embedded sources are 
found. The curves include both the efficiency of iden- 
tifying the embedded galaxies with SExtractor and the 
correlation and cross correlation selection cuts. Below 
0.4 " galaxies are merged by SExtractor and the pair is 
lost. For dim galaxies, which tend to have higher red- 
shifts, noise in the galaxy detection becomes an increas- 



Simulations are used to estimate both the number of 
lensed galaxy pairs based on the local density of galaxies 
in the COSMOS survey as well as the efficiency of find- 
ing those pairs. The simulation of lensed galaxy pairs is 
the same as used previously in our GOODS search [l5| . 
The idea is to use the density of galaxies in our catalog 
of sources to monte carlo the number of pairs that would 
exist from any theoretical string crossing our survey. An 
important aspect of this calculation is the simulation of 
the galaxy redshifts. We use the same parameterization 
of the redshift distribution as [25] . For these studies, 
COSMOS tile 55 is used as a typical example. The sim- 
ulated redshift distribution of the galaxies in this tile is 
shown in Fig. [3l 

The advantage of a catalog-level signal simulation is 
that we can quickly embed as many strings as needed 
to get an accurate estimate of the average number of 
lensed galaxies observed from a particular set of string 
parameters, including redshift and £sin/3. An example 
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FIG. 4: Efficiency of detecting pairs of galaxies lensed by a 
cosmic string as a function of pair opening angle and redshift. 



5 




String Redshift G^/c 2 (10" 3 ) 



FIG. 5: 95% upper limits for lensed galaxies produced by a 
cosmic string as a function of the string mass and redshift. 
Dashed lines are the individual limits from each angular bin 
shown above. The solid line is the average limit for all direc- 
tions. 



ingly important effect. 



IV. RESULTS: 

The distribution of matched galaxy pairs shown in 
Fig. [2] rises nearly linearly. This is expected for the back- 
ground pairs and any excess above the background could 
be evidence of a cosmic string. For comparison, simulated 
pairs from a cosmic string at a redshift of 0.5 and S sin j3 
of 4 " are included in the top panel. They are normalized 
to the mean length of a string crossing the survey, 1.19°. 
The upper curve is the total simulated signal. The lower 
curve includes the detection inefficiencies from Fig. 3J 

We compare a wide variety of predicted cosmic string 
signals to the data to determine limits. Signal pairs, 
n s , are summed from 0.4 " to a maximum opening an- 
gle beyond which there is no more signal. The observed 
matched pairs, n G 5 S , and the background, n^g, from 
Fig. [2] are summed over the same range of opening an- 
gles. We then compute the single-sided Neyman 95% 
confidence limit, n^ m , which is the minimum signal that 
is consistent with background fluctuations. Signals with 
n s > are excluded by the data. The resulting 95% 
upper limits are shown in Fig. [5] All four angular bins 
yield similar limits that extend from 1 " < Jsin/3 < 7". 
We average the four angular bins for the global limit. 
Taking the mean tilt of a string with respect to the ob- 
server to be < sin/3 >= 2/tt we relate the opening angle 
to the mass scale via the factor 87r^ = S < sin/3 > 
shown on the right-hand axis. We see no evidence for 
cosmic strings out to a redshift greater than 0.5 and place 
a 95% upper limit of G/i/c 2 < 3.0 x 10~ 7 . 



FIG. 6: 95% upper limits on the global parameter, Q s trv, 
as a function of the string mass. 



If strings are rare occurrences, it is likely that none 
would appear in the COSMOS field (fov = 1.57 degrees 2 ) 
and that other surveys may yield different results. Be- 
cause strings are line-like objects, we expect that the 
COSMOS field of view would intersect a single string 
about 0.5% of the time. Using a simple geometric monte- 
carlo we determine that approximately 820 straight 
strings randomly placed in a volume of comoving radius, 
77, are needed for a 95% detection rate in a COSMOS 
fov. In this case, the average COSMOS fov contains 3 de- 
tectable string crossings as expected from Poisson statis- 
tics. This string density corresponds to a total length 
of string, L to t = 1200?], in the volume and we use the 
invariance of L tot /r] to set limits on fl strings- 

_ K L tot/y)y 8ttG 

strings ~ (4/3)^ 3 3tf 2 [ ] 

where 77 is the comoving distance computed with h = 0.7, 
tt M = 0.27, and ^ A = 0.73. Fig. [6] shows the string 
densities excluded by this method. The limit excludes a 
string density that is 0.28% of the critical density for the 
smallest mass strings and rises to 0.4% for more massive 
strings. Each limit assumes that all cosmic strings have 
the same Gji/c 2 as predicted by some models. Because 
non-uniform string networks are predicted to have a very 
steep spectral index, (-10<index<-6) [5], strings at the 
measurement threshold will dominate the sensitivity. 

A wide variety of theoretical models predict T = 
^strings/ (87T Gji) of order 10 [5[. Our measurement is con- 
sistent with T = 375 for low-mass strings and V = 110 for 
our largest-mass strings and thus indicates that we are 
not yet sensitive to masses of current theoretical interest. 
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V. CONCLUSION: 

We have used the COSMOS survey field to search for 
cosmic strings. We find no evidence for the gravitational 
lensing signature. We have included the observational 
efficiencies in our analysis using the same technique we 
used previously on the GOODS survey. Figs. [5] and 
[6] summarize our results. We set 95% upper limits of 
Gfi/c 2 < 3.0 x 10 -7 out to redshifts greater than 0.6 
which leads in turn to a global limit of ^strings < 0.0028. 
We note that, while these results have as their foundation 
the hypothesis of long straight strings, they also exclude 
strings with moderate curvature. 

The global limit on ^strings is nearly 10 times stronger 
than our previously published limit [15]. We want to em- 
phasize that our technique is complementary to other 
methods. We excluded masses that are smaller than 
those excluded by other direct CMB searches [11[ . How- 
ever, our masses are larger than those reported by pa- 



rameter fits to the CMB @4Iq| as well as gravitational 
wave searches [12|, [l3|- The recent microlensing search 
carried out on SDSS quasars [16| is sensitive to very low- 
mass strings, but reports a similar sensitivity for Q strings- 
Our direct searches may be sensitive to a particular class 
of strings that could be missed by other searches. 
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